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Abstract 

In our version of Watts and Strogatz's small world model, space is a (i-dimensional 
torus in which each individual has in addition exactly one long-range neighbor chosen 
at random from the grid. This modification is natural if one thinks of a town where 
an individual's interactions at school, at work, or in social situations introduce long- 
range connections. However, this change dramatically alters the behavior of the contact 
process, producing two phase transitions. We establish this by relating the small world 
to an infinite "big world" graph where the contact process behavior is similar to the 
contact process on a tree. 

Keywords: contact process; small-world network; phase transition; epidemic 

1 Introduction 

Small world graphs were first introduced by Watts and Strogatz (1998). In their model, 
they take a one-dimensional ring lattice and connect all pairs of vertices that are distance 
m or less. They then "rewire" each edge with probability p by moving one of the ends 
at random, where the new end is chosen uniformly. This leads to a graph that has small 
diameter but, in contrast to the Erdos-Renyi model, has a nontrivial density of triangles. 
These are both properties that they observed in the collaboration graph of film actors, the 
power grid, and the neural network of the nematode, C. elegans. 

The small world model has been extensively studied, although most investigators have 
found it more convenient to study the Newman and Watts (1999) version in which short- 
range connections are not removed and each neighbor has a connection to a long-range 
neighbor with probability p (See Figure 1). Its graph theoretic properties (e.g., the average 
distance between two points and the clustering coefficient) are well understood (see Albert 
and Barabasi (2002) for the physicist's view point or Barbour and Reinert (2001) for rigorous 
results). Our focus here will be on the behavior of processes taking place on these networks. 
Chapter six of Watts(1999) discusses the SIR (susceptible-infected-removed) disease model 
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on small world graphs in which individuals that are susceptible (state 0) become infected 
(state 1) at a rate proportional to the number of infected neighbors. Infected individuals, 
after a random amount of time of fixed distribution, become removed (state 2), i.e., forever 
immune to further infection. 

The SIR model on the small world graph has a detailed theory due to its connection to 
percolation: we draw an oriented edge from x to y if x will succeed in infecting y during 
the time it is infected and the persistence of the epidemic is equivalent to percolation. See 
Section 8.2 of Newman (2003) for what is known about the SIR models on small world 
graphs. Here we will investigate the more difficult SIS (susceptible- infected-susceptible) 
model, known to probabilists as the contact process, where recovered individuals are im- 
mediately susceptible. Moreno, Pastor- Satorr as and Vespigiani (2002) have studied this 
model by simulation, but we know of no rigorous results for the contact process on the 
small world. Berger, Borgs, Chayes, and Saberi (2004) have proved rigorous results for the 
contact process on the preferential attachment graph of Barabasi and Albert (1999). 

Our version of the small world will be as follows. We start with a d-dimensional torus 
(Z mod R) d and connect all pairs of vertices within distance m of each other using the 
|| • | |oo norm. We require R to be even so that we can partition the R d vertices into R d /2 
pairs. Consider all such partitions and then pick one at random. A new edge is then drawn 
between each pair of vertices in the chosen partition. When m <C R, we think of these new 
edges as long-range connections. We will call this graph S^, keeping in mind that this is a 
random graph. 

The reason for insisting that all individuals have exactly one long-range neighbor is that 
we can define an associated "big world" graph B m that is non-random. We define to 
consist of all vectors ±(21, . . . , z n ) with n > 1 components with Zj G Z d and Zj 7^ for 
j < n. Neighbors in the positive half-space are defined as follows: a point +(z\, . . . , z n ) 
is adjacent to +(z\, ... ,z n + y) for all y with < ||y||oo < m (these are the short-range 
neighbors of +(z\, . . . , z n )). The long-range neighbor is 



See Figure 2 for a picture of the one dimensional case with m = 1. Of course in this case 
the graph is a tree. 

The drawing of the small world in Figure 2 is more convenient for our proof, but the 
graph can be more succinctly described if all the long-range edges point down. In this case 
it is the free product Z d * {0, 1} where the second factor is Z mod 2. Elements of this group 
have the form zq1z\\ ...\z n where Z{ G Z d — {0} for < i < n. In words, this is the point 
you reach by moving zq in the first copy, going down a long-range edge, moving sideways 
by zi, going down a long-range edge, etc. 
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We will consider the discrete-time contact process. On either the small world or the 
big world, an infected individual lives for one unit of time. During its infected period it 
will infect some of its neighbors. All infection events are independent, and each site that 
receives at least one infection is occupied with an infected individual at the next time. A 
site infects itself or its short-range neighbors, each with probability a /(2m + l) d . It infects 
its long-range neighbor with probability (3. Let A = a + (3 and r = a/ (3. Hereafter we will 
assume a > (3, and we fix the ratio 1 < r < oo. We will use Bt to denote the contact 
process on the big world and £t for the contact process on the small world. 

The number of sites within distance n of a given site in the big world grows exponentially 
fast, so it is natural to guess that its contact process will behave like the contact process on 
a tree. Consider a tree in which each vertex has the same degree, let be a distinguished 
vertex (the origin) of the tree, and let A$ be the set of infected sites at time t on the tree 
starting from occupied. For the contact process on the tree or on the big world, we can 
define two critical values: 

Ai = inf{A : P(|A?| = eventually) < 1} (1) 
A 2 = inf{A : liminf P(0 6 A?) > 0}. 

t^oo 

We call Ai the weak survival critical value and A2 the strong survival critical value. Pemantle 
(1992) showed that for homogeneous trees where every vertex has at least four neighbors, 
Ai < A2. He and Liggett (1996) who extended the result to trees with degree 3, did this by 
finding numerical bounds on the two critical values which showed they were different. Later 
Stacey (1996) found a proof that did not rely on numerical bounds. For another approach 
to separating the critical values see Lalley and Selke (2002). 

Theorem 1.1. For each < r < 00 there exists too such that for all m > too, Ai < A2 for 
the contact process on B m . 

The proof of Theorem 1 1 . 1 1 will follow from Propositions ll.ll and ll.2l below. Our approach 
to proving Theorem II .11 is to find an upper bound U (Proposition II . l)l on the limiting weak 
survival critical value lim sup^^^ Ai , a lower bound L (Proposition ll.2|) on the limiting 
strong survival critical value liminfm^oo A2, and then to show that U < L. 

To obtain the two bounds, we will compare the contact process to the branching random 
walk which has births like the contact process, but allows more than one particle at a site. 
More explicitly, each individual lives for one unit of time and sends offspring to its site and 
to its short-range neighbors with probability a/(2m + l) d and to long-range neighbors with 
probability (3. Birth events are independent, and the state of a given site on the next time 
step is the total number of births there. 

The number of particles of the branching random walk at a given site stochastically 
dominates the number of particles of the contact process at a given site. Therefore weak 
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(strong) survival of the contact process on B m implies weak (strong) survival of the branch- 
ing random walk on B m . However, in order to get an upper bound on the weak survival 
critical value of the contact process, we need the opposite implication. This is achieved by 
showing that in the limit as m — > oo, the behavior of the contact process is essentially the 
same as that of the branching random walk. The following proposition, proved in Section |5J 
shows that the trivial necessary condition for weak survival gives the asymptotic boundary 
of the survival regime. 

Proposition 1.1. If a + (3 > 1 then the contact process survives for large m. 

To obtain a lower bound L on A2, we use the fact that strong survival of the contact 
process on B m implies strong survival of the branching random walk on B m . Let \^ w (m) 
be the strong survival critical value of the branching random walk. To compute the limit 
of \ b 2 W {m), we define the "comb" of degree m, C m , by restricting B m to vertices of the 
form {+(z), +(z, 0), — (0)} and all edges between any of these vertices. As before, +(z) 
and +(z,0) are long-range neighbors as are +(0) and —(0). The short-range neighbors of 
+(z) are +(z + y) for < ||y||oo < m - The vertices +(z, 0) and —(0) have no short-range 
neighbors. To see the reason for the name look at Figure 3 which gives a picture of the 
graph for m = 1, d = 1. As we will show in Section |31 linim^oo A?,"" (to) is the critical value 
for survival of the branching random walk on the comb. 

It will be shown in Section |3] that the branching random walk on the comb survives if 
the largest eigenvalue of 




is larger than 1 (since the comb is like Z rf the weak survival and strong survival critical 
values are the same). Solving the quadratic equation (a — A)(— A) — f3 2 = the largest root 
is 

a + ^Ja 2 - 4/3 2 
2 

A little algebra shows that this is larger than 1 exactly when a + (3 2 > 1 . The proof of the 
following proposition is given in Section 03 Since the critical value for the branching random 
walk is smaller on the comb than it is on the big world, this result is sharp asymptotically. 

Proposition 1.2. If a + fi 2 < 1 then there is no strong survival in the contact process for 
large m. 

Comparing the above with Proposition 11.11 shows that we have separated the critical values 
for all (3 < 1. 

We have not been able to generalize Stacey's elegant argument to the big world graph. 
However, we have been able to establish the following ingredient used in his argument, 
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which was proved for the tree by Morrow, Schinazi, and Zhang (1994). Here, and in what 
follows, we will write as shorthand for +(0). 

Theorem 1.2. There exist constants C\ > and C2 such that 

exp(C 2 t)<E(|£ t °|)<-Lexp(C 2 t). 

Moreover, C2 < when A < Ai, C2 = when A = Ai, and C2 > when A > Ai. 

Let re = min{t : B® = 0} be the extinction time of the contact process on the big world. It 
follows from Theorem II . 21 that rg < 00 with probability one when A < Ai and when A < Ai, 
rg has a geometric tail. 

Having established the existence of two phase transitions on the big world, our next 
question is: How does this translate into behavior of the contact process on the small 
world? Let ag = min{i : B® = or 6 B®} be the first time that the infection either 
dies out or comes back to the origin starting from one infection there at time 0. When 
Ai < A < A2, t~b and as are both infinite with positive probability, and when A > A2, o~b is 
almost-surely finite. Let t$ = min{i : $ = 0} and as = min{t > 1 : £° = or € be 
the corresponding times for the contact process on the small world. 

Theorem 1.3. Writing => for convergence in distribution as R — ► 00 we have 
(a) ts is stochastically bounded above by tq and t$ => rg. 
(6) as is stochastically bounded above by a& and as =^ a&. 

Theorems 11.11 and 11.31 show that the small world contact process exhibits different be- 
havior in the three regimes (0, Ai), (Ai,A2), and (A2,oo). In the first interval the process 
dies out. In the second, a& = 00 with positive probability, so the time for the infection to 
return to is not tight (in the sense of weak convergence). In the third interval, erg < 00 
with probability one, so the infection returns to in a time that has a limit as the system 
size increases. This is especially interesting since other researchers have studied the contact 
process on the small world by simulation without having noticed this qualitative difference 
between (Ai,A2), and (A2,oo). 

Since the small world is a finite graph, the infection will eventually die out. However, 
by analogy with results for the (i-dimensional contact process on a finite set, we expect that 
if the process does not become extinct quickly, it will survive for a long time. Durrett and 
Liu (1988) showed that the supercritical contact process on [0, R) survives for an amount 
of time of order exp(ci?) starting from all ones, while Mountford (1999) showed that the 
supercritical contact process on [0,R) d survives for an amount of time of order exp(cR d ). 
At the moment we are only able to prove the last conclusion for the following modification 
of the small world contact process: each infected site infects its short-range neighbors with 
probability a/(2m+l) d and its long-range neighbor with probability (5, but now in addition, 
it infects a random neighbor (uniformly) from the grid with probability 7 > 0. 
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From a modelling point of view, this mechanism is reasonable. In addition to long-range 
connections with friends at school or work, one has random encounters with people one sits 
next to on airplanes or meets while shopping in stores. In the language of physics, the 
model with 7 = has a quenched (i.e., fixed) random environment, while the model with 
(3 = has an annealed environment. 

Our strategy for establishing prolonged survival is to show that if the number of infected 
sites drops below rjR d , it will with high probability rise to 2r]R d before dying out. To do 
this we use the random connections to spread the particles out so that they can grow 
independently. Ideally we would use the long-range connections (instead of the random 
connections) to achieve this; however, we have to deal with unlikely but annoying scenarios 
such as all infected individuals being long-range neighbors of sites that are respectively 
short-range neighbors of each other. 

Theorem 1.4. Consider the modified small world model on a d-dimensional torus of radius 
R with random infections at rate 7 > 0. If X > Ai and we start with all infected individuals 
then there is a constant c > so that the probability the infection persists to time exp(ci? d ) 
tends to 1 as R — > 00. 

This result is somewhat surprising and is qualitatively different from the one Stacey (2001) 
proved for a finite homogeneous tree with radius R. In that case if one starts with all sites 
occupied, the extinction time grows linearly in R for A < A2, but grows doubly exponential 
in R for A > A2. In contrast Theorem 11.41 is. to our knowledge, the first result on any 
graph for which the survival time of the contact process is exponential in volume when 
Ai < A < A2. Of course the reason we can get such a result is because a truncated tree is 
much different from our locally tree-like small world. 

We should mention that Berger, Borgs, Chayes, and Saberi (2004) have studied the 
contact process on Barabasi and Albert's (1999) scale-free graph, which has a power law 
degree distribution ~ Ck~ 3 and a critical value A c = 0. For all fixed A > 0, they were able 
to show survival up to time exp(cA 2 n 1//2 ) for the contact process on a graph with n vertices, 
but again, this is sub-exponential in volume. To do this they combined an estimate for the 
survival time of the contact process on a star graph with large degree, with the fact that 
the largest degree in their random graph is 0(n 1 / 2 ). 

Acknowledgments. The authors would like to thank a referee for suggesting that 
the results be formulated in terms of the two parameters a and 0, an associate editor 
for indicating the connection between the big world and the free product, and Laurent 
Saloff-Coste for a tutorial on free products. 
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2 An Upper Bound on Weak Survival 



To prove Proposition 1.1, consider the branching random walk Bi on the half-space 13^ 
consisting of all points of the form +(z±, . . . , z n ). We say a particle is at level / if it is at a 
point of the form +(z±, . . . , z\). 

Clearly, the branching random walk is supercritical when a + (3 > 1. The idea is to 
run the supercritical process to time t large enough so that we get enough special particles 
at different points all of the form +(zi, . . . , 0) on some level I. We can then consider 
the half-spaces formed by all points at level greater than or equal to / treating points 
+(zi, ... ,zi— i,Q) as the origins. Each of the special particles starts its own branching 
random walk on these disjoint half-spaces, and by construction none of these branching 
random walks will ever run into each other (although a given branching random walk can 
run into itself). If the expected number of special particles at time t is greater than 1, then 
a comparison with a branching process shows that the total number of special particles on 
levels nI,iiGNat times nt is infinite with positive probability. We then choose m large so 
that the contact process looks like a branching random walk. 

Let us now carry out this plan more rigorously. Letting 

( 1 \ 
a = a [ 1 — — -t-7 , 

we have that E(|B t |) > a{a + f3) l ~ l . One way to see this is by disallowing short-range births 
onto points of the form, +(z\, . . . , z^m-ii 0). The particles in B t are on levels 0, . . . ,t so 
there is a level £(t) where the expected number of particles is at least a(a + + 1). 

When a + (3 > 1 we can choose m and T large enough so that 

a{a + l3) T - 1 _3_ 
T+ 1 > &P' 

Modify the dynamics so that the only births occurring at time T + 1 are those across 
short-range edges on level £{T) (excluding births to points of the form +{z\, . . . , z^m-i; 0)), 
and the only births occurring at time T + 2 are those across long-range neighbors from level 
£(T) to level l(T) + 1. These birth at time T + 2 will lead to special particles on level 
£(T) + 1, and the half-spaces obtained by eliminating the long-range edges just crossed are 
all disjoint. 

The expected number of special particles in the modified branching random walk on 
level £{T) + 1 at time T + 2 is at least 3. We claim that if m is large then the expected 
number of special particles on level £(T) + 1 at time T + 2 in the contact process on the 
half space is at least 2. 

The expected number of particles in the branching process at any time t < T is at most 
(a + 0) T . Let E K be the event that \B t \ < K(a + [3) T for all t < T. For e > choose K 
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large enough so that 

E(|B T |l £jr )>(l-e)E(|Br|). 

We now compare the contact process to the branching random walk. To go from time t 
to time t + 1 in the contact process, we first let births occur across long-range edges. This 
cannot lead to two generation t + 1 particles landing on the same site. We then let the 
short-range births occur one at a time. On the event Ex, the probability of a collision on 
any given birth is at most K(a + (3) T /(2m + l) d , so the probability of a collision on some 
birth at some time t < T is at most TK 2 (a + [3) 2T /(2m + l) d . 

Now taking m large enough, the expected number of special particles at level £(T) + 1 at 
time T + 2 in the contact process is at least 2. Since these particles sit atop half spaces that 
are disjoint, the number of special particles on level 2(i(T) + 1) at time 2(T + 2) for these 
processes in their disjoint half-spaces are independent. By comparing with a branching 
process we conclude that the contact process survives. 



3 A Lower Bound on Strong Survival 

The first step is to compute the critical value for the branching random walk on a comb. 
This is not necessary for the proof of the lower bound on strong survival but it explains the 
expression for the bound that we obtain. Additionally, it motivates us for the latter part 
of this section. Here we will return to our original notation and fix r = a//3 > 1 and vary 
\ = a + p. 

Lemma 3.1. The critical value for the branching random walk on the comb C m is indepen- 
dent of m and is equal to 

2(r + l) 
r + Vr 2 + 4 

Proof. We consider sites +(z) to be level 1 and all other sites on the comb to be level 2. 
The level transition probabilities for a random walk on the comb are given by 

/ _a_ \ / _r_ _j_ \ 

I o+/3 a+13 _ I r+1 r+1 ) 

[a o ° ) 

Solving the quadratic equation 

— x I (— x) — —77 = x 2 -x — — tt , (2) 

r + 1 ) K ' (r + 1) 2 r + 1 (r + 1) 2 w 

the largest eigenvalue of the matrix is 

r + \A" 2 + 4 
2(r + l) 
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This shows that the expected number of particles at tends to exponentially fast if 

» - 2(r+1) 

r + yr 2 + 4 

and gives the lower bound on the critical value. For the other direction, we note that if the 
random walk is on level 1 after n steps then there have been m down steps, m up steps, and 
n — 2m sideways steps. The sideways steps are independent of the up and down steps so 
the probability of being at after n steps, given that the walk is on level is greater than 
cjn d l 2 . This shows that the expected number of particles at grows exponentially fast if 

w 2(r + l) 



r + \Jr 2 + 4 

The proof can now be completed using techniques of Madras and Schinazi (1992) or Pe- 
mantle and Stacey (2001). □ 

Since the branching random walk r)f on B m dominates the branching random walk t]f 
on the comb of degree m, we have 

X b 2 rW (C m ) > \ b 2 rW (B m ) 

However, this inequality is in the opposite direction of what we want. To get a result in the 
other direction, we introduce another graph structure Km which is similar to B m except 
that has no short-range neighbors, and we replace copies of Z rf by copies of the complete 
graph on M = (2m + l) d vertices. Again, each point can be described algebraically by a 
vector (#!,..., z n ) with n > 1 integer components, but now we have z\ = and < Zj < M 
where Zj ^ for 1 < j < n. Figure 4 shows a picture of Km with M = 5. 

Consider two branching random walks r]f and 77^ on the respective graphs B m and Km 
starting from one particle at the origin. In the second process, individuals that would be sent 
to short-range neighbors of stay at 0. We can couple rjf and rfc so that rjf(0) < r)t(0) 
for all t. 

Let \2 TW {Km) denote the strong survival critical value of rjf '. In order to get a lower 
bound L for the strong survival critical value in terms of A and r, it is enough to show that 

2(r + 1) . . . . j ru , 



< liminf X^ W (K M )- (3) 

If we then interpret this in terms of a and /?, Proposition II. 21 will have been proved. 

Let Sk be the random walk on K m which jumps to short-range neighbors with probability 
a/ ((2m + l) d (a + 0)) and Ion g-range neighbors with probability f3/(a + (3). Jumps from 
to short-range neighbors make no change in the location of the particle. Define 

r = inf{£; > 1 : S k = 0} and F x = ^ A fe P(r = k). 

k>l 
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Lemma 3.2. X b { w {K, M ) > sup{A : F x < 1}. 

Proof. For the branching random walk started with one particle at 0: 

E(r?f (0)) = X k P(S k = 0). 

Summing over k we define 

G A = £A fc P(S fc = 0). 

fc>0 

We need only show that F A < 1 implies G A < oo, since this in turn implies that the 
branching random walk dies out locally. Breaking things down according to the value of r 
we have for k > 1 

k 

p(s fc = o) = £p(T = OP(s fc _, = o). 
i=i 

Multiplying by X k and summing over 1 < k < oo we get 

G x - 1 = F A G A 

so that G A = 1/(1 - F x ). □ 

To examine the behavior of S k we will compare it to a birth and death chain = <f>(Sk) 
on the nonnegative integers, (j) maps in KLm to 0, while for n > 2, a point of the form 
(zi, . . . , z n _i, 0) maps to 2n — 3 and a point of the form (z±, . . . , z n ) with z n / maps to 
2n — 2. See Figure 4 for a picture of the mapping. Let u = 1/(1 + r) and M = (2m + l) d . It 
is clear from the symmetries of /C m that i?fc is a Markov chain with transition probabilities: 
r(0, 0) = 1 — u and r(0, 1) = u. When j is odd, the last coordinate is so: 

r(j,j-l) = u, r(j,j) = (l-u)/M, r(j, j + 1) = (1 - «)(1 - 1/M). 
When j > 2 is even the last coordinate is non-zero so: 

r(j,j + l)=u, r(j,j-l) = (l-u)/M, r{j, j) = [l - u){\ - l/M). 
It follows from our definitions that r = inf{& > 1 : R k = 0}. For x > 0, let 

/ l (x) = ^A fc P a; (r = A:) 

fc>i 

and define /i(0) = 1. By considering what happens on one step, if x > 

h(x) = \J2r(x,y)h(y). (4) 

y 
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Using this equation for x = 2n, 2n — 1 and 2n + 1, we have for n > 1 

fc(2n)^l - = >?H 2 n - 2)r(2n, 2n - l)r(2n - 1, 2n - 2) 

+ A 2 /i(2n + 2)r(2n, 2n + l)r(2n + 1, 2n + 2) 

+ h(2n) Xr(2n, 2n) ^1 - A ^ ~ ) + A 2 r(2n, 2n - l)r(2n - 1, 2n) 
+ A 2 r(2n,2n+ l)r(2n + l,2n) . 
This can be rewritten as ah(2n + 2) — bh(2n) + ch(2n — 2) = where 
a = A 2 «(l - «)(1 - 

„ _ 1 _*fl^_ (A _*'(! -) )(1 -„)(!- i,) 

_ A 2 (l-t*) 2 (l-£) _ a ^ 2 
M 

X 2 u{l - u) 
° ~ M ' 

The solutions to the homogeneous difference equation are of the form h(2n) = C\Q™ + 6*2^2 
where 6\ and 62 are the roots of a6 2 — b6 + c = 0, or 



b + Vb 2 - Aac , . 6 - V^ 2 - 4ac 

6*i = and 9 2 = - • 

la la 

When A < 1, h(2n) is decreasing so h(2n) = 9%- Since h(x) is an analytic function of A, 
this is the formula for all A inside the radius of convergence. When M — > 00, c — > while 
a and b have positive limits so using the Maclaurin expansion of y/1 — x we get 



e 2 = b Y 

2a b 

and it follows that h{2) — » 0. 

Intuitively what we have shown is that particles that reach level 2 (i.e., fall off the comb) 
and their descendants can be ignored. To complete the calculation now we observe that 
using @ with x = 1 and letting M — » 00 gives h(l) — > Xu. By considering what happens 
on the first step starting from we have 

F x = A(l - u + uh(l)) -> A(l - u) + u 2 A 2 . 
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Setting F x = 1 and recalling u = 1/(1 + r) gives the quadratic equation 



' A 2 + A — - — -1 = 0. 



(1+r) 2 r + 1 

The change of variable A = 1/x gives (J2J and establishes © • 

4 Exponential Growth and Decay 

In this section we will prove Theorem 11.21 Since we have been writing for +(0), we will 
similarly write —0 for —(0), and —1 for one of the nearest neighbors of —(0). We begin 
with a well-known property of trees, see e.g., Lemma 6.2 in Pemantle (1992). 

Lemma 4.1. Let T 3 be the rooted binary tree in which the root has degree one and all other 
vertices have degree 3. Suppose there are I infected sites on T 3 . There must be at least 1 + 2 
copies o/T 3 ; disjoint except for possibly the root, where the only infected site is the root. 

Proof. We proceed by induction. If / = 1 there are 3 = 1 + 2 copies. Each time we add a 
vertex we destroy one copy and create two more. □ 

See Figure 5 for a picture of a case with 1 = 9. We have drawn B\ but it is the same 
as T 3 . As predicted there are 11 copies of T 3 , marked by *'s on the edges leading to their 
roots. We will call these edges, "exterior edges." The tree result extends to the following 
property of the big world. 

Lemma 4.2. Given a finite set B C B m , there is a set of vertices D C B with \D\ > ^-\B\ 
such that each vertex x G D is either (a) adjacent to a vacant translate of a copy of B^ 
with x sitting at — 0, or (b) adjacent to a vacant translate of a copy of £>+ U {—0} with x 
sitting at a site that is a nearest neighbor of— 0. 

Proof. Consider first the case d = 1. We can embed T 3 into B m . Using Lemma [4. H it follows 
that there are at least (I + 2)/3 completely disjoint copies of T 3 where the only infected 
site is the root. If the exterior edge from the boundary vertex x is a long-range edge then 
we are in case (a). If the exterior edge from the boundary vertex x is a short-range edge 
we are in case (b). See Figure 5. 

In d = 2 we can still embed a tree but it has a variable degree. In each plane we connect 
(0, 0) to its four nearest neighbors. For k > 1 we connect (k, 0) to both (k + 1, 0) and (k, 1) 
while we connect (k, j) to only (k, j + 1) for j > 1. Then extend the construction to the 
other four quadrants so that it is symmetric under 90 degree rotations. See Figure 6 for a 
picture. All points in addition have long distance neighbors so the result holds with C2 = 5. 
We leave the details for d > 3 to the reader. □ 
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With Lemma f4.2l established. the rest of the proof is very similar to the proof of Theorem 
2 in Madras, Schinazi, and Zhang (1994). 

Lemma 4.3. For the contact process B® on B m , there exists C\ > such that 

E(|<J) > CiE(|B°|)E(|B,°|). 

Proof. For j = 0, 1, let \Bf' 3 \ be the number of points in B^ for the process started with 
a single infected individual at —j. Let Mt = rain{~E(\B^~'~°\), E(|S^"' 1 |)}. By additivity, 
the Markov property, and Lemma 14.21 we have 

E(|<J) > -E(\B?\)M S . 

Cd 

The lemma will be proved when we show that there exists C% > such that M s > C\~Ei(\B®\). 
To do this we begin by observing that for j = 0, 1 

where we are considering the event that on the first step — j infects —0, and on the second 
step —0 infects itself and 0. The factor of 1/2 comes from the fact that when starting from 
{—0,0}, the expected number of points in the positive half-space is one-half the number in 
the whole space. To take care of the differences in starting configuration and time, we note 

E(\B°\) <E(\B^- 0} \) < \ 2 B(\Bl%- 0} \). 

From this it follows that for j = 0, 1 

liminf 5M3>><?V <• 



t^oo E(|5 t °|) ~2 \(2m + l) d X / 

□ 

Proof of Theorem M.'A By the additivity of the contact process and the transitivity of B m , 

E(|S° +S |) < E(|£ t |)E(|.B°|). 
Subadditivity then implies that 

lim ilogE(|fl t °|) = inf ilogE(|B t °|) = C 2 (6) 

t^oo t t>0 t 

exists so that exp(C2i) < E(|i? t °|). A similar argument using superadditivity of CiE(|S°|) 
and Lemma 14.31 shows that 

lim ilogCiE(|£?°|) = sup^logdEd^l) = C 2 (7) 
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so that ^-exp(C 2 t) > E(|B t °|). 

We turn now to the sign of C2 in the different regimes. When A > Ai we use the 
well-known fact that on the event 

{|.B t | > 1 for all k > 0}, 

lim^oo \ Bt\ = 00 almost surely. So by the upper bound ^- exp(C 2 /c) > E(|B t °|), we see that 
C 2 > when A > Ai. 

When A < Ai, the work of Aizenman and Barsky (1987) gives 

£E(|£ t °|)<oo 
t>o 

so that C2 < in this regime. 

For A = Ai we use a continuity argument. Note that for any fixed t > 0, 

A-ilogE(|itf|) 

is a continuous function of A > 0. By (jSJ) and (J7|), A — ► C2 is lower and upper semi-continuous 
and thus continuous. Therefore when A = Ai it must be that C2 = 0. □ 

5 Proofs of Theorem 11.31 and 11.41 

Proof of Theorem M.ft We start by viewing B m as a "covering space" of 5^. Given a 
realization of the random graph <S^, we mark a distinguished vertex which we identify 
with the origin of B m . From there we identify the long-range edge of the distinguished 
vertex with the long-range edge of the origin in B m . Similarly we identify the short-range 
edges of the two graphs together in such a way so that graph distances on the underlying 
structure given by B m are preserved. Continuing in this manner we can identify each vertex 
in B m with some vertex in our realization of S^. 

Using this identification, couple together the processes B® and £® for each realization 
of the random graph S^. By additivity of the contact process, it is clear that \B® | > \^\ 
which shows that t$ is bounded above by t&. Choose e > and N > 0. Using the graph 
identification above, let Gr^n be the event (on the random graph S^) that there is no cycle 
of length 2N which contains 0. For fixed N, we have that P(Gr jv) — ► 1 as R — > 00 (if 
this is not clear to the reader, an explicit argument is given in the proof of 11.41 below). If 
there are no cycles of length 2N containing in the graph then no two vertices within 
distance N of are identified together in B m . 

Choose i?o so that R > Ro implies P(G/j i at) > 1 — e. Using the coupling to identify 
points of B m with S^, we have for all R > Rq, 

P(5£ = Ck for all k < N) > 1 - e 
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which proves lim^oo ts = Tg. The proof of limjj_ J . o0 as = o~b is similar and is therefore 
omitted. □ 

Proof of Theorem \1.4\ Let M = (2m + l) d be the number of neighbors and let 

8M 

The contact process on the big world with 7 = is supercritical so we can pick T so that 
> 5/5 and we can pick K so that if £y is the contact process restricted to the ball 
of radius K around then E(|^|) > 4/5. Let Nk be the number of points in the ball of 
radius K around on the big world. Pick rj so that 

77 ■ < — 

' 4M ~ 5 

It suffices to show that if the number of occupied sites drops below r]R d occupied sites 
then with probability greater than 1 — ex.p(—bR d ) for some b > 0, it will return to having 
more than 2r\R d occupied sites. To see this we just choose c < b giving us a high probability 
of having exp(cn d ) successful recoveries. Since each recovery takes at least one unit of time, 
the result follows. 

The first step towards achieving the above is to show that the first time it drops below 
r]R d it does not fall too far. To get a lower bound on how far it falls, we only look at 
births from sites onto themselves. These are independent events with probability a/M. 
A standard large deviations result for the binomial implies that if there are at least r\R d 
particles alive at time t then the probability of less than r\R d a/2M alive at time t + 1 is 
less than exp(— b\R ) for some 61 > 0. Thus the probability this occurs at some time before 
exp(6ii? d /2) is less than exp(-6ii2 d /2). 

After the number of occupied sites falls below r]R d , our next step is to randomize the 
locations by having one time step in which we only allow births induced by the parameter 
7. If we start with r]R d a/2M particles, then after these particles give birth onto their ran- 
domly chosen neighbors (with no other births allowed) we will have, on average, rjR d ja/2M 
particles. In fact, if t + 1 is the time at which the number of occupied sites falls below r\R d 
(but not too far), then using the large deviations estimate for the binomial once more we 
have that the probability of there being less than r\R d ^a / 4M particles at time t + 2 is less 
than exp(— b2R d ) for some 62 > (we will take care of double counting below). Having 
used the randomized births in this step, we ignore them for the rest of the proof and use 
the process with 7 = which we have assumed is supercritical. 

We say that a site x is good if the small world is identical to the big world inside a ball of 
radius K around x. Fixing x, the probability of a self-intersection (i.e., the probability that 
x is not a good site) when we generate the ball of radius K around x is less than N^/n. 
To see this, we grow the ball around x starting from just x and its short-range neighbors. 
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Adding the long-range neighbor of x, we see that the probability it causes a self-intersection 
is bounded above by Nx/n. But we also have to be worried about the short-range neighbors 
of this long-range neighbor causing a self-intersection, so we increase the bound generously 
to N^/n. Now since there are at most Nk neighbors to add, the bound N^/n follows; 
therefore, the probability that x is a bad site tends to as n — > oo. 

Now fix a realization of the small world graph such that the fraction of bad points does 
not exceed 1 /20. We want to find a subset G of the randomized births that lie on good sites 
and for which the evolutions of £f ,x £ G are independent. To do this we imagine that the 
successful randomized births occur one at a time (successful with respect to 7). We accept 
the first birth if it lands on a good site. We accept the second birth if it lands on a good 
site and it does not fall within the ball of radius 2K around the first birth, and we continue 
in this manner. 

When there are j successful randomized births accepted, the total number of sites that 
lie within the balls of radius 2K around them does not exceed jN2K- When 

j < 7]R d ~fa/4M, 

the probability of the next successful randomized birth landing in the forbidden zone is less 
than 1/5 by our choice of r\ above. Since we are on a realization of whose fraction of bad 
points does not exceed 1/20, the probability of discarding a successful randomized birth 
conditioning on j prior acceptances is less than 1/4. Thus the number of births discarded 
out of the first r}R d ^a/ (8m + 4) successful randomized births is bounded by a binomial with 
success probability p = 1/4. Our large deviations estimate implies that with probability 
less than exp(— bsR d ), the number of randomized births remaining, \G\, is larger than 

8M 1 

By our choice of G, the particles in G can grow independently up to a distance of K, 
or in other words, the processes £f,x G G evolve independently. Since E(|£y,|) > 4/(5, a 
final large large deviations estimate for bounded i.i.d. random variables tells us that the 
probability of ending up with fewer than 2rjR d particles at time T is less than exp(— b±R d ). 
The total of our error probabilities is less than exp(— bR d ) which completes the proof. □ 
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Figure 1. The Newman and Watts version of a small world graph. For simplicity we have 
not drawn the short-range connections. 
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Figure 2. The big world graph, B\. 
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Figure 3. The comb C±. We have rotated the edge from +(0) to —(0) to make the picture 
better match the name. 
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Figure 4. Comparison graph /Cm when M = 5. Numbers next to the vertices indicate the 
corresponding states in the birth and death chain. 
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Figure 5. Finding exterior edges (*'s) on the big world graph, B. Black dots indicate 

infected sites. 
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Figure 6. Embedding a tree in Z 2 . 



24 



